Interqubit coupling mediated by a high-excitation-energy quantum object 
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We consider a system composed of two qubits and a high-excitation-energy quantum object used to 
mediate coupling between the qubits. We treat the entire system quantum mechanically and analyze 
the properties of the eigenvalues and eigenstates of the total Hamiltonian. After reproducing well- 
known results concerning the leading term in the mediated coupling, we obtain an expression for 
the residual coupling between the qubits in the off state. We also analyze the entanglement between 
the three objects, i.e. the two qubits and the coupler, in the eigenstates of the total Hamiltonian. 
Although we focus on the application of our results to the recently realized parametric-coupling 
scheme with two qubits, we also discuss extensions of our results to harmonic-oscillator couplers, 
couplers that are near resonance with the qubits and multi-qubit systems. In particular, we find that 
certain errors that are absent for a two-qubit system arise when dealing with multi-qubit systems. 



INTRODUCTION 



Superconducting qubits are among the main can- 
didates for the possible implementation of quantum- 
information-processing tasks Coherent dynamics of 
a single qubit has been achieved at various laborato- 
ries. Several interestingtwo-qubit experiments have also 
been performed 0, 0,1, i, 0, 0, i, ©, El, El El, El, 




Wi . Il9l [2fj|. The early experiments 
ilillflllil El were limited to fixed 
inter-qubit coupling. In order to scale up qubit circuits, 
however, it is highly desirable to be able to tune in situ 
the coupling strengths between the different qubits. The 
idea of coupling two qubits to a high-excitation-energy 
quantum object that would mediate coupling between 
the qubits [H, [H, [U HE HI has lead to experi- 
mental demonstrations of tunable coupling [H, EH EH • 
However, Refs. [H, E3 probed the magnetic properties 
of the circuit in its ground state; thus the approach 
used there is not suited to implement the standard gate- 
based quantum computing, but possibly adiabatic quan- 
tum computing [27| . where the quantum register is ide- 
ally never excited. In Ref. [l5[ tunability of the coupling 
was demonstrated through spectroscopic measurements 
with the qubits biased away from their coherence opti- 
mal points. Generalizing the idea of mediated coupling 
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|25| . |26| ] to parametric tunable coupling 
it was proposed that one can bias the 
qubits at their optimal points and also adjust the dc 
component of the mediated coupling to cancel the direct 
inter-qubit coupling throughout the experiment. Apply- 
ing a microwave pulse to the coupler at the sum or dif- 



ference frequency of the qubits' characteristic frequencies 
would then turn on the coupling, but only for the dura- 
tion of the applied microwave pulse. This proposal, com- 
bining tunable coupling and long coherence times, was 
realized experimentally in Ref. [171 ]. 

Another related direction of growing research activity 
is the idea of using a harmonic-oscillator 'cavity' as a data 
bus with the potential that a single cavity could mediate 
coupling between a large number of qubits [l9l I20I l32l. |33| . 
If the harmonic oscillator mediates coupling via high- 
energy virtual excitations, it can be described similarly 
to other high-excitation-energy couplers. Couplers that 
are near resonance with the qubits, but sufficiently de- 
tuned such that they mediate coupling through virtual 
excitations, can also be treated using a similar approach. 
We shall show, however, that it would be rather diffi- 
cult to achieve tunability in the coupling in the case of 
harmonic-oscillator couplers. 

Early theoretical studies on tunable couplers have gen- 
erally taken the semiclassical approach, which we shall 
explain below. The semiclassical treatment is sufficient 
to evaluate the leading term in the effective coupling me- 
diated by the coupler. As qubit circuits that include 
couplers are becoming an experimental reality, however, 
there is an increasing need for a more careful analysis 
of these quantum couplers. Some recent studies [29L l34j 
used a number of quantum-mechanical techniques and 
obtained results beyond the semiclassical calculations. 
For example, Ref. [3J] derived an expression for the resid- 
ual coupling term when the main coupling channel is 
turned off. Here we take a fully quantum approach, 
where we analyze the properties of the eigenvalues and 
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eigenstates of the total Hamiltonian. We obtain results 
that were not captured by previous studies, particularly 
results concerning the ideality of the off state in this cou- 
pling scheme. 

The main questions of interest that the quantum treat- 
ment can be used to answer include: (1) Can we make the 
coupler's energy splitting very large but still maintain the 
mediated coupling between the qubits? (2) Can we turn 
the coupling off completely? In other words, will there 
be any residual coupling terms if the system is biased 
such that the main coupling channel is at its zero point? 
(3) How much entanglement is there between the qubits 
and the coupler, and how does this entanglement affect 
a realistic experimental setup? These questions will be 
answered in the analysis below. 

The paper is organized as follows. In Sec. II we de- 
scribe the system and its Hamiltonian. In Sec. Ill we 
review the semiclassical approach and discuss what pre- 
dictions we can expect from it regarding the ideality of 
the off state. In Sec. IV we perform the fully quantum 
analysis of the problem: We derive expressions for the 
leading-order effective coupling strength mediated by the 
coupler, the residual coupling in the off state and the 
amount of entanglement between the qubits and the cou- 
pler; We also discuss the extension of our results to the 
case of a harmonic-oscillator coupler and that of a coupler 
that is almost resonant with the qubits. We discuss in 
Sec. V the implications of the results obtained in Sec. IV 
with regard to present-day and future experiments. In 
this context we also consider multi-qubit systems. Sec- 
tion VI contains concluding remarks. Some details of the 
calculations are explained in the Appendix. 



II. SYSTEM AND HAMILTONIAN 

Let us take a system composed of two qubits and a 
third object that we would like to use as a coupler. Al- 
though tunability is considered the main advantage of 
this coupling scheme, for the purpose of answering the 
questions of main interest to us it suffices to focus on the 
case where the external bias parameters are set to fixed 
values. We therefore treat a time- independent Hamilto- 
nian. Since we are assuming that we can speak of three 
distinct quantum objects, we must be able to write down 
a Hamiltonian that reflects this clear separation of the 
different objects in the system. We therefore express the 
Hamiltonian as: 



h 2 = f^(2) + fa^(2) 

2 2 



H = H i + Hi + He + H12 + Hic + H. 



2Ci 



(1) 



where the first three terms are the Hamiltonians of the 
separate objects in the system, and the last three terms 
describe coupling between those objects. As a realistic, 
representative case, we take the different terms in the 
Hamiltonian to have the forms: 
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with the coupler energies 0, rji, r]2, ■■ 
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and similarly for B. For the case of superconducting 
flux qubits, Aj is the minimum gap of qubit j, ej rep- 
resents the deviation from the degeneracy point of half- 
integer flux quantum threading the qubit loop [i.e. ej = 
/ pj ($ ex tj — $o/2), where I p j is the persistent current 
of qubit j, & cx tj is the externally applied flux in the 
loop of qubit j, and $ is the flux quantum], J is the di- 

(7) 

rect qubit-qubit coupling strength, and <Tq ,' are the usual 
Pauli matrices of qubit j. Note that the minimum gap 
Aj is the coefficient of oiP above, in contrast with some 
alternative conventions used in the literature where Aj is 

the coefficient of ai^ . Note also that the operators A and 
B must satisfy the relations Aij — A*,- L and Bij — B*; L . 

The coupler's Hamiltonian He is written in its own eigen- 
basis (thus it is diagonal), and its ground-state energy 
has been set to zero. We shall express Hamiltonians and 
energies in frequency units throughout this paper. 

We now make the assumption that the largest energy 
scale in the problem is the excitation energy of the cou- 
pler. In other words r\i is larger than any relevant energy 
scale in the Hamiltonian excluding He- We shall not 
make any assumption regarding the relation between the 

qubit energy scale \J A 2 j + e 2 and the qubit-couplcr cou- 
pling energy scale (i.e. the scale of the matrix elements 
in A and B) . Although it is not crucial for most of our 
analysis below, we shall generally take the direct coupling 
energy scale J to be smaller than the qubit energy scale. 



III. SEMICLASSICAL TREATMENT 

Before embarking on the semiclassical description of 
our system, it is instructive to recall the problem of cal- 
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culating interatomic forces within the hydrogen molecule 
f35l ]. One starts by assuming that the nuclei have fixed 
locations in space, leaving the degrees of freedom asso- 
ciated with electron motion as the only variables in the 
problem. The ground state of the electronic degrees of 
freedom is calculated and expressed as a function of the 
distance between the nuclei. At this point, the total en- 
ergy (direct Coulomb energy of the nuclei plus electronic 
ground-state energy) is expressed as a function of the 
relative position operator between the nuclei, and the 
nuclear motion is treated quantum mechanically. One 
can now obtain information related to atomic motion 
or molecular states without having to worry about elec- 
tronic motion. The first step in the calculation ensures 
that the effects of the electrons' adiabatic adjustment 
to nuclear motion are properly taken into account. The 
reason why the electronic degrees of freedom can adjust 
adiabatically to nuclear motion is that they are associ- 
ated with a much higher energy scale, and thus changes 
that result from nuclear motion are felt by the electrons 
as very slow variations. 

A similar procedure can be applied when dealing with 

inter-qubit couplers. Since the qubit-coupler interactions 

(i) 

contain the operators <xir , one first assumes that the 
qubits are in eigenstates of dx . One therefore needs to 
calculate the ground state energy of the Hamiltonian (or, 
more precisely, the four Hamiltonians): 

Hc,eS = He + H\c + H 2 c 

er (i) =±l - 

H C ±A±B. (9) 

The ground state energy [i.e. the four values obtained 
from Eq. ©] can then be expressed in the form (36[: 

£0 {#C,eff} = ci + c 2 a { ^ + c 3 <jW + aa^ ® <4 2) . (10) 

Note that we have started with the assumption that the 

(7) 

qubits have well-defined values of tri . One could there- 
fore say that the above expression is not a quantum op- 
erator. However, it is straightforward to follow the adi- 
abaticity argument: the coupler will, to a very good ap- 
proximation, always be in the ground state that corre- 
sponds to the instantaneous values of a4 (including the 
possibility of quantum superpositions). One can there- 
fore turn back to the qubits and analyze their dynamics 
using the effective Hamiltonian: 

-ffq,cff = Hi + H 2 + Hl2 + -^mediated, (H) 



where 

^mediated = Ci + CaffW + C 3 4 2) + <**P ® ^ ■ (12) 

The above derivation shows that the response of the cou- 
pler to external perturbations (induced by the qubits and 
represented by the terms ±A and ±-B) translates into 
a renormalization of the qubit bias points [second and 
third term in Eq. (|12[) ] and an additional coupling term 
between the qubits [last term in Eq. ([12])]- 

It is interesting to note here that in the semiclassical 
treatment above, the effect of the coupler on the qubits 
is completely described by the Hamiltonian i? mo diated- 
Therefore, if the bias point is chosen such that the coef- 
ficient C4 in Eq. (TT2]) cancels the direct coupling strength 
Jo, there would be no residual coupling between the 
qubits. The off state would therefore correspond to com- 
plete decoupling between the qubits. In contrast, we shall 
show below that the fully quantum treatment predicts 
the presence of finite residual coupling effects. 



IV. QUANTUM TREATMENT 

We now take the Hamiltonian of Sec. II (Eq. [T]) and 
treat it quantum mechanically. In particular, we would 
like to analyze the properties of the eigenvalues and 
eigenstates of the total Hamiltonian. For the purposes 
of the calculations in this section, we divide the Hamil- 
tonian into two parts as follows: 

H = H Q + V, (13) 

where 

Ho = ^ + ^+H G (14) 

V = + y4 2) + H 12 + flic + H 2C . (15) 

We use the basis {|000), |010), |100), 1 1 10) , |001), |011), 
1 101) 1 1 111) , 1 002), ...}, where the first, second and third 
quantum numbers describe, respectively, the state of the 
first qubit, second qubit and coupler (the state of the cou- 
pler will be distinguished from those of the qubits using 
an underline throughout this paper); note that we allow 
the coupler to have more than two levels, and we shall 
assume an infinite number of levels in the expressions 
below. We can now express Hq as 
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A. Coupling strength: leading term 

In order to calculate the effective coupling strength, 
we now calculate the energies of the lowest four energy 
levels. We therefore want to construct a 4 x 4 effective- 
Hamiltonian matrix describing the lowest energy levels 
while taking into account the effects of the higher levels. 



We now follow a standard calculation [33] (see Ap- 
pendix) that gives: 
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(20) 



where the operator P projects the state onto the space of 
the lowest four eigenstates of H (|000), |010), 1 100) and 
1 110)); alternatively, one could say that the operator P 
removes the size mismatch between the four-dimensional 
Hilbert space of interest and the infinite-dimensional op- 
erator V . The sum over I in Eq. (|20|) runs over states 
where the coupler is in one of its excited states. The en- 
ergies Ek,i are the eigenvalues of Hq with the qubits in 
state k and the coupler in state I. In order to find the 
lowest-order expression for the mediated coupling term, 
we make the approximation H e g — Ek,i ~ — f]i in Eq. (|20|) . 



Thus we obtain the expression 
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where 



A s hift = 



1=1 



\A 0l \ 2 + \B 0l \ 



(22) 



(2) 

The overall shift i? shift does not have any physical con- 
sequences to this order of the calculation and can be ne- 
glected. 

The effective Hamiltonian can now be expressed as 



where 



H c ft — H o + ^coupling, 



H Q = ^M + |^) + ^(2) + £^(2) 



(23) 



2 " z ' 2 x 2 

-'-'coupling — u u x ^ u x ' 



2 * 



with the parameters 



ei = ei + 2A 00 
£2 = £2 + 2£?oo 
J = Jo-Jl. 



(24) 



The above results agree with those of Ref. [29( when the 
parameters of our model are taken to correspond to those 
considered in Ref. [29| . 

First we mention the physical interpretation of the 
terms 2Aqq and 2Bqo in the expressions for ej. Taking 
the experimentally relevant case of flux qubits as an ex- 
ample, these terms describe the fluxes generated by the 
coupler (in its ground state) and going through the qubit 
loops. As a result, if one wishes to bias the qubits at their 
optimal points, the externally applied fluxes through the 
qubit loops are not set to $o/2, but they are shifted from 
that value to compensate for the coupler-induced contri- 
butions. This is the physical explanation of the difference 
between €j and ij [381 ]. 

We can now answer the question of how high the cou- 
pler's energy splitting can be. We note that the mediated 
coupling strength J\ is second order in the scale of A and 
B, and it is inversely proportional to the scale of r/i. Al- 
though this means that we cannot simply take r\\ — > oo 
keeping the other parameters fixed (otherwise J\ — * 0), it 
also means that if we increase A nm , B nm and r\i (keeping 
the ratio A iBi /r]i fixed) we can in fact take 771 — > 00 
while maintaining the same level of mediated coupling. 
We shall see below that this situation is desirable for 
purposes of reducing residual-coupling effects. 

We are also in a position to comment on the question 
of monostability of the coupler (i.e. the idea that a single 
energy level of the coupler is relevant in the system un- 
der consideration) . The presence of the matrix elements 
Aiq and Biq, which describe coupling between the cou- 
pler's ground and excited states, in the expression for the 
mediated coupling strength (Eq. |2"2"1) demonstrates that 



the excited states of the coupler play an important role 
in the coupling mechanism. The semiclassical treatment 
relies on the fact that the relevant information contained 
in these matrix elements is also encoded, and more easily 
accessible experimentally, in the response of the coupler's 
ground state to weak perturbations. As a result, knowl- 
edge of the matrix elements themselves is not necessary; 
knowledge of simple response parameters is sufficient in 
order to calculate the mediated coupling strength. 



B. Residual coupling in the off state 

In Sec. IV. A we have derived the leading-order terms 
in the effect of the coupler on the two-qubit system. The 
results agree with the predictions of the semiclassical 
treatment (although no explicit expressions were given in 
Sec. III). Our main motivation for using the fully quan- 
tum treatment, however, is that it allows us to go further 
in the calculation, e.g. evaluating any off-state residual 
coupling terms in the two-qubit effective Hamiltonian. 
Such closer examination of the off state will be the sub- 
ject of this and the following subsection (as well as part 
of Sec. V). 

Since we are interested in an ideal off state with the 
qubits at their optimal points, we consider the case where 
J — e\ = £2 = (we shall allow J to take nonzero val- 
ues in the course of the discussion below). In order to 
enhance the robustness of the decoupling between the 
qubits, we take Ai 7^ A2. The effective Hamiltonian 
(Eq. [23]) is diagonal to lowest order in this case. We can 
therefore proceed with calculating higher-order correc- 
tions using the states |000), |010), |100) and 1 1 10) . 

The calculation of the residual-coupling Hamiltonian is 
now carried out by calculating the energies of the lowest 
four levels of the entire system (we shall refer to them as 
£-000, -Eflio, -E10Q1 and Eno). The four values that we ob- 
tain for the energies in the four-level spectrum can then 
be used to extract four quantities: (1) an overall energy 
that can be neglected, (2,3) the corrected (i.e. renormal- 
ized) values of the qubit splittings Aj and (4) a residual- 
coupling energy (given by E 00 q + -E110 - ^010 - £100). It 
is this last quantity that is of most interest to us here. It 
characterizes a coupling Hamiltonian of the form: 



-^residual — ^residual & 



(1) 



where 



Jvi 



En 



En 



E, 



010 



E 



100 



ddual 



(25) 



(26) 



Note that in writing Eq. (|25p we have used our assump- 
tion that the qubits are biased at their optimal points. 

Although one can make analytic progress calculating 
the energies using the pseudo-degenerate perturbation 
theory approach, reaching the relevant results requires 
several steps of the iterative procedures explained in 
the Appendix (the lowest-order corrections to the in- 
dividual energies do not contribute to the combination 
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-Eboo + £110 — E010 — E100). We therefore only present 
the results of numerical calculations that find the energy 
levels of the entire system. It also turns out that proceed- 
ing with the rather general model used so far complicates 
the extraction of the important results. We therefore 
have to make some simplifying assumptions. From now 
on, we focus on a simple case that also happens to be 
relevant to recent experiments: We take the coupler to 
be a two-level system, such that 



Hr 



H2C 



£25.(0) + ^c^(c) 
2 2 
Ac 



A c / 

— - [tan 6c 

Jic ^ 
J 2C ® 40) 



J(C) 
(C) 



r (C) 



(note that, unlike the convention used in Sec. II, 
not diagonal here). We therefore have 771 = \J A c - 



(27) 

Hq is 

,2 _ 



A c /|cos6» c |, A m = -An = -Jicsin^c, A w = A m = 
Jic cos^c and similarly for B. As a result, we find from 
Eq. 



that 



J = Jo 



2JicJ2C cos 3 6 C 



(28) 



For clarity and definitcness in the following analy- 
sis, we now focus on the parametric-coupling scheme 
[2a . I29I [30| . and we start by noting three factors that 
can contribute to determining the ideal dc bias point (it 
is worth mentioning at this point that in general it will 
not be possible, nor necessary, to satisfy all three con- 
ditions simultaneously). Firstly, we note from Eq. (|28]) 
that the range of values for the effective coupling strength 



J extends from J — 2JicJ2c/Ac (at 6c = 0) to J 
(at dc — tt/2). In order to maximize the effective 
coupling strength of the parametric-coupling scheme, it 
would be desirable to set the dc bias point such that 
J = Jo — Jic»/2c/Acj i-e. halfway between the two 
extremes; This situation would be obtained by setting 
6c ~ 0.65 ~ 7r/5, irrespective of the specific system pa- 
rameters. At this bias point, the allowed amplitude of the 
ac driving signal (i.e. before encountering nonlinearities 
in J) is maximized (see Fig. la). Secondly, it is desir- 
able to set J = at the dc bias point. Thirdly, it is also 
desirable to set the residual coupling energy J res iduai = 0. 

We now proceed with the numerical calculations as 
follows: We first fix the parameters Ai, A2, Ac and Jo 
(e.g. Ai = 4 GHz, A 2 = 5 GHz, A c = 30 GHz and 
Jo = 50 MHz). The coupling strengths Jic and J 2 c 
are then chosen as Jic = J 2 c an d JicJ2c/Ac = Jo- 
This choice ensures that a single dc bias point satisfies 
at least two of the three desirable conditions mentioned 
above, i.e. J vanishes while maximizing the achievable 
effective coupling strength for the parametric-coupling 
scheme. The coupler's bias point, determined by the an- 
gle 6c, is now treated as the only variable in the problem. 
For each value of dc , we (numerically) bias the qubits at 
their optimal points (determined by simultaneously min- 
imizing both qubit energy splittings) and calculate the 
residual coupling energy. The results are shown in Fig. 1. 
We can clearly see that there is some residual coupling 
in the off state (i.e. at 6c ~ 0.65); however, the residual 
coupling strength in Fig. 1 is very small compared to typ- 
ical relaxation and dephasing times. After varying all the 
parameters by at least a factor of 2 in either direction, 
we can identify that the residual-coupling term is given 
by 



7 A1A2J0J1CJ2C Sq AiA 2 Ji C J| c 5 . 2 , . , . 

/residual ~ 4 ^ cos D 6 C - a U cos 6 C sin (20c), (29) 



where the coefficient a is approximately 20. 

It is encouraging for future experimental efforts that 
the residual-coupling term obtained above decreases 
rapidly with increasing coupler gap Ac, even if the me- 
diated coupling strength is kept at a fixed level. Further- 
more, it is possible to adjust the system parameters such 
that the conditions of vanishing J and J rC siduai are simul- 
taneously satisfied. An example of this situation is shown 
in Fig. 2. Note that the achievable effective coupling 
strength for the parametric-coupling scheme has been 
reduced from that in Fig. 1, because the dc bias point 
has now moved away from the desirable point 6c ~ 0.65 
mentioned above. 

Another relevant quantity when analyzing the ideal- 
ity of the off state is how the effective coupling strength 



changes in response to changes in the coupler's bias point. 
This quantity can be derived easily from Eq. (|28| : 

dJ _ 6 Jic J2C cos 4 6 C sin 6 C , , 

dec' A c • m 

With proper design parameters, this quantity also de- 
creases rapidly with increasing Ac- 

C. Entanglement in the energy eigenstates 

In addition to the residual-coupling Hamiltonian, an- 
other natural question to ask when analyzing the ideality 
of the off state is how much entanglement there is be- 
tween the qubits and the coupler, and between the qubits 
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6 c 

FIG. 1: Coefficients of the main (a) and residual (b) coupling 
terms (i.e. J and Jresiduai) as functions of the coupler's bias 
point 8c- In generating this figure we took Ai =4 GHz, 
A 2 = 5 GHz, A c = 30 GHz, J = 50 MHz and J 1C = J 2 c = 
V Jo Ac- It should be noted, however, that the results are 
only weakly dependent on the exact choice of parameters. 

themselves, in the different energy eigenstates. Unlike 
the residual-coupling Hamiltonian, where we had to re- 
sort to numerical calculations, this entanglement can be 
evaluated using the analytic approach. The analysis of 
how the energy-eigenstate entanglement manifests itself 
as errors in a quantum calculation will be postponed until 
Sec. V. 

Here we are interested in the nearly ideal situation 
where any residual entanglement is treated as a small 
error. We therefore take the situation analyzed in the 
previous subsection ( J = ?i = e 2 = 0, Ai ^ A 2 , and the 
coupler is a two-level system) , and we analyze the degree 
of entanglement to lowest non- vanishing orders. 

As explained earlier in this section, we can start our 
calculation using the states |000), |010), 1 100) and 1 1 10) . 



For the lowest order of this calculation, we can use the 
pseudo-degenerate perturbation theory calculation ex- 
plained in the Appendix. We find that 



IV'oo) ? 


a |000) - 


2r]i 


Mo 
Vi 


lioi) 


-^|011) 
Vi 


I^Ol) * 


a |010) - 


^lioo)- 

2r]i 


Mo 

Vi 


1111) 


_ |001) 

m 


l^io) p 


a |100> - 


^loio)- 

2r?i 


Mo 

Vi 


1001) 


-^°|ni) 
m 


iV'n) * 


* |no)- 


A |000) - 
2/71 ' 


Mo 

Vi 


ion) 


-— 1101X31) 

Vi 



up to a normalization constant slightly smaller than 1. 
Note that the above expressions can be collectively sum- 
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FIG. 2: Same as in Fig. 1, except that Jic = Jic = 
0.75V Jo Ac- As opposed to Fig. 1, J and Jresiduai now vanish 
at the same bias point. 



marized as (39j : 



\ipnm) ~ \nm0) —\nmQ) —\nml) —\nml) 

1 _ g, 42) _ ^10 .a) & (0 _ ^(2) 3 ai c)\ |nmQ) 

2T]i 1]i 7)1 J 

1 - J -^W a ^i 2) - *£2±°°W ) - ^!f£,( 2 ) j |ronfl>> (32) 



Two notes are in order here: (1) The mixing of computational-basis states in the energy eigenstates 
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(Eq. [32)) cannot be identified as arising from a simple 
residual-coupling Hamiltonian. In other words, one can- 
not find a two-qubit Hamiltonian that reproduces any of 
the corrections in Eq. (|3"2")h This can be perhaps most 
clearly seen in the fact that the signs in front of the 
small terms in Eq. (f32|) are independent of the state. 
In contrast, one always expects different states to ac- 
quire corrections of different signs when dealing with a 
direct coupling Hamiltonian. Higher-order effects in this 
system therefore cannot be completely described by sim- 
ply using small corrections to the reduced (i.e. 4x4) 
two-qubit Hamiltonian. (2) The above expressions were 
obtained assuming only that the coupling strength J van- 
ishes, ?i = ?2 = and Ai ^ A2. They are not affected by 
the value of the residual-coupling Hamiltonian discussed 
in Sec. IV. B. This fact demonstrates that the vanishing 
of both J and Jresiduai does not imply that the qubits are 
entirely decoupled. Corrections to the energy eigenstates 
still have to be considered when calculating possible er- 
rors, as will be explained in Sec. V. 

As in the case of the residual-coupling Hamiltonian, it 
is encouraging for scalability considerations that the en- 
tanglement in the energy eigenstates decreases and ap- 
proaches zero with increasing Ac, assuming that Jo and 
Ji are kept fixed. One can therefore say that the cou- 
pler becomes more and more ideal as Ac is increased. It 
should be noted that the superconducting gap for alu- 
minum is ~ 50 GHz. It would be interesting to investi- 
gate in the future how serious a constraint this number 
imposes on the maximum allowed value of Ac- 



D. Harmonic oscillator as a coupler 

We can use the results derived above for a two-level 
coupler to infer the corresponding results for a harmonic- 
oscillator coupler. We now consider the situation where 



term 



He 








H2C 





(33) 



where a) and a are the creation and annihilation oper- 
ators for the coupler. The above coupling Hamiltonians 
change the state of the coupler by one excitation. It 
should therefore be a good approximation to truncate 
the coupler to its lowest two energy levels. We now find 
that 



J= Jo 



2 Jic J2C 



and 



J rcsid 



ual 



AiA 2 Jo Jig 'he 



(34) 



(35) 



In order to study the tunability of a harmonic- 
oscillator coupler, we now consider the effect of an ap- 
plied external field on the mediated coupling, adding the 



H 



field 



h(a + a f ) 



(36) 



to the Hamiltonian. By defining b = a + 2h/uj(j, we can 
see that the new Hamiltonian takes the same form it had 
in the absence of Hfi^, except that a and aJ are replaced 
by b and b' , and the qubit bias parameters ei and £2 are 
shifted to ei and e 2 (just as discussed in the general case 
above). In particular, the values of J and J res iduai are 
not affected by the applied field. The mediated coupling 
is therefore not tunable, assuming that the coupler's bias 
point is set by an applied linear field, i.e. a field that 
affects the coupler according to Eq. (|3"6"|) . 

If the frequency of the harmonic oscillator (or alterna- 
tively Jic and J2C) were tunable, one would be able to 
obtain a tunable value of J (see Eq. IMl) . Possible de- 
signs for such tunable couplers have been proposed the- 
oretically [4dl | . but they have not been realized experi- 
mentally. Given the advantages they could provide in 
terms of tunable coupling, it would be highly desirable 
to fabricate such tunable oscillators in the future. 

A related design would be to use a (possibly non- 
tunable) harmonic-oscillator element, which can be made 
large in size, in the circuit and add a tunable coupler be- 
tween this oscillator and each qubit around it. With this 
design one would gain the advantages of both (1) the 
large 'cavity' being able to function as a data bus con- 
necting a large number of qubits and (2) the coupling 
being tunable. With this design one would also be able 
to use the parametric-coupling scheme to perform entan- 
gling operations on any qubit in the circuit and the cavity 
relatively easily (using red- or blue-sideband transitions), 
with the qubits biased at their optimal points. 



E. A coupler that is near resonance with the qubits 

So far, we have focused on the case where the cou- 
pler's excitation energy is the largest energy scale in the 
problem. For completeness we consider in this subsection 
the case where the coupler's excitation energy is close to 
those of the two qubits, but with sufficient detuning to 
avoid exciting the coupler. In other words, we consider a 
situation similar to the experiment of Ref. [20l ]. but not 
that of Ref. 

The largest energy scale in the problem is now the 
qubit and coupler energy splittings, which are almost 
equal. The next largest energy scale is the detuning 
between the qubits and the coupler (this assumption 
is made in order to avoid large entanglement between 
the qubits and the coupler in the energy eigenstates, in 
which case we would have to deal with real excitations 
in the coupler). The detuning between the qubits, direct 
interqubit coupling strength and qubit-coupler coupling 
strength can take any values. 

Because the qubit energy scale is the largest energy 
scale in the problem, the mediated-coupling term now 
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describes an excitation transfer between the two qubits. 
As such, the form of the effective-coupling term will only 
make sense in connection with the single-qubit Hamilto- 
nians, not the physically defined operators a x and a z as 
in the case of a high-excitation-encrgy coupler. Assum- 
ing that the qubits are biased at their optimal points and 
the coupler is a harmonic oscillator with only one rele- 
vant excited state (and using the notation of Sec. IV. D), 
the mediated-coupling term for this excitation-transfer 
process is well described by the Hamiltonian 

fj J\cJ2C /_(1)_(2) , \ 
-^mediated = X \ °+ °- + h.C. I , (37) 

ujc — A V / 

where A = (Ai + A2V2 (note that |Ai — A2I is assumed 
to be much smaller than \ojc — A|; note also that the 
qubit splittings Aj will be slightly modified because of 
the interaction with the coupler). Using numerical cal- 
culations we find that the residual-coupling term is given 
by Eq. fl25|) with 

T _ JoJichc ,„o\ 

^ residual . - \ 2 ' V / 

(wc - A) 

The fact that the residual-coupling strength scales in- 
versely with the second power of the qubit-coupler de- 
tuning as opposed to the fourth power of the coupler's 
energy splitting indicates that residual coupling can be a 
more serious issue in this case than in the case of a high- 
excitation-energy coupler. It should be noted, however, 
that in the experiment of Ref. [2(| Jo was essentially zero, 
in which case residual coupling should be negligible. 

V. ESTIMATING ERRORS IN A TYPICAL 
EXPERIMENT 

In this section we discuss a typical present-day ex- 
perimental procedure. Using concrete examples, we will 
be able to discuss rather clearly the possible errors that 
the residual coupling and energy-eigenstate entanglement 
might cause. 

In most, if not all, present-day experiments, the qubits 
are controlled using a single microwave line. This situ- 
ation, however, cannot be maintained for larger multi- 
qubit systems, where the energy levels of the entire sys- 
tem become densely packed. We shall therefore assume 
that local control lines are used to address the individual 
qubits. 

A. Two-qubit system 

The usual theoretical approach to describing a quan- 
tum system that is composed of several distinct objects 
commonly runs along the following line of reasoning: 
First we assume that these objects are in a separable 
initial state, then we controllably evolve the system us- 
ing the total Hamiltonian, we take the trace over the 



degrees of freedom of the unused objects (in this case the 
coupler), and we find the answer to the physical ques- 
tions of interest. Given the entangled form of the energy 
eigenstates obtained in Sec. IV, one might expect that 
this 'contamination' of the states will reduce the observed 
coherence effects. Indeed, following the above-described 
recipe, one finds such a reduction in observable coherence 
effects, e.g. reduced gate fidelities. As we now explain, 
however, this is not the correct description of a typical 
experiment, and the resulting errors and limitations are 
also described incorrectly in this way. 

First let us consider the assumption about the ini- 
tial conditions. In a real two-qubit experiment (as- 
suming essentially zero temperature), the system starts 
in its ground state \tpoo) 1 000). When preparing 
the desired initial state (e.g. a product of two single- 
qubit states), single-qubit operations are performed us- 
ing small-amplitude pulses that are resonant with indi- 
vidual qubits. These pulses are typically weak enough 
that the single-qubit operations are performed over times 
that are long compared to l/|Ai — A2I (the energy dif- 
ference |Ai — A2I is used here because it is, to a good 
approximation, the smallest energy difference in the spec- 
trum of the four lowest energy levels). Rather than 
perform ideal single-qubit operations, such weak pulses 
drive transitions between the eigenstates of the Hamil- 
tonian of the entire system, regardless of the form of 
these eigenstates. One could therefore ignore the fact 
that \ipnm) 7^ |nm0) and simply use the eigenstates of 
the Hamiltonian (i.e. IV'oo), \1p01), IV'io) and |i/>n)) as 
the computational basis for performing a quantum al- 
gorithm. One can then use weak pulses with properly 
calibrated frequencies to perform any operation on this 
effective two-qubit system. Thus the initial and all sub- 
sequent states in the experiment can be very close to the 
desired form, except that the basis states do not have the 
simple, separable form [e.g., the state (IV'oo) + \^w))/V^ 
is prepared instead of the state (|000) + |100))/%/2]. 

There are two main possible sources of errors that can 
affect the above picture. They are both related to the 
fact that in most quantum algorithms one needs to per- 
form operations on a certain qubit without knowing the 
states of the other qubits. The residual-coupling Hamil- 
tonian (Eq. |25|) provides the first source of errors. For 
example, the resonance frequencies for the transitions 
IV'oo) IV'io) an d IV'oi) IV'n) are shifted from each 
other by a frequency difference given by 4J res iduai- As 
a result, the resonance frequency for performing opera- 
tions on qubit 1 depends on the state of qubit 2, and 
vice versa. In principle, one could perform a single-qubit 
gate using two pulses, with each pulse being resonant 
with one of the transitions. This approach, however, 
is impractical for many-qubit systems, where resonance 
lines generally split into an exponentially large number 
of lines. A more scalable alternative is to perform single- 
qubit operations using driving amplitudes that are large 
compared to the spread in the relevant resonance fre- 
quencies, such that the pulse can be considered on res- 
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FIG. 3: The error-estimating quantity <557R a bi,i = 

((V'oo|a"i 1) |i/'io} — (i/>oi\ffx \ipu))/2 f° r tne relative spread in 
the Rabi frequency of qubit 1 as a function of the coupler's 
bias point 9c for the same parameters as in Fig. 2. Similar 
results would be obtained for qubit 2. 

onance for all the relevant transitions. However, the 
residual-coupling Hamiltonian will then cause an unde- 
sirable phase accumulation over the course of running the 
algorithm. Until this residual coupling is suppressed in 
future experiments (e.g. by increasing Ac as discussed 
in Sec. IV. B), one might need to apply refocussing pulses 
to reduce its effects. As can be seen from Figs. 1 and 
2, pushing J ros iduai to the kHz range should be possi- 
ble using realistic experimental parameters. The other 
possible source of errors is the fact that the Rabi fre- 
quency of qubit 1 oscillations generally depends on the 
state of qubit 2, and vice versa (This splitting of the Rabi 
frequencies is a result of the mixing in the energy eigen- 
states) . Such errors can be described by quantities of the 

form 5fi Ra bi,i = ((^io|o-£ 1 i^oo) - (i>n\vx \ipoi))/2. It is 
straightforward to see that, to the lowest order given in 
Sec. IV.C, 

(^ol^l^oo) = (Vn|4 1) IV'oi). (39) 

In order to obtain an estimate for <5f2R a bi,ij we perform 
some numerical calculations. We show in Fig. 3 the re- 
sults for the parameters of Fig. 2. With these parameters 
we can see that ^Rabu ~ 10~ 8 . This type of errors can 
therefore be made extremely small using realistic exper- 
imental parameters. 

A third, but less serious, point of potential concern re- 
lates to quantities of the form (V'oo |<5"4 1 ' 1 IV'oi}- This quan- 
tity describes, in some sense, how much qubit 2 'feels' 
a signal applied to qubit 1 (Similarly to the splitting of 
Rabi frequencies mentioned above, this type of error is a 
result of the mixing in the energy eigenstates). Using the 

parameters of Fig. 2, the error estimator (ipoo\o'i 1 ' > \ipoi) 
is of order 10~ 2 , and vanishes close to the point where 
J = 0. This type of errors can be suppressed further by 
using microwave amplitudes such that the small fraction 
of a signal applied to qubit 1 that is felt by qubit 2 is 



small enough to be considered far off resonance. Alter- 
natively, one could say that single-qubit operations on 
qubit 1 must be performed on a time scale that is long 
compared to | (V'oo 1 0x |?/>oi)|/|Ai — A 2 |. This is a very 
mild constraint, even for existing experiments. 

Finally, a tricky issue related to the usual theoreti- 
cal approach mentioned at the beginning of this subsec- 
tion is when the coupler degrees of freedom are traced 
out at the end of the actual quantum calculation. Here 
again, the measurement process in typical experiments 
does not follow the simple picture of a sudden, accurate 
measurement device that probes the basis 0/1. For flux 
qubits, for example, one typically designs the readout 
device to simply be a very accurate device for measur- 
ing magnetic fields (measuring the flux generated by the 
qubits and going through the loop of a readout SQUID). 
If the qubit-SQUID coupling is weak compared to the 
qubits' energy scale and it acts for a sufficiently long du- 
ration, the readout device will perform the measurement 
in the basis of the eigenstates of the system Hamiltonian, 
not the physically defined clockwise and counterclockwise 
current states. If, in addition, the relevant decoherence 
rates in the SQUID are small enough, the readout fidelity 
will not be limited by the exact amount of mixing be- 
tween the clockwise and counterclockwise current states 
in the energy eigenstates [4l[. Therefore, the readout 
device can, in principle, give the correct reading (0 or 1) 
essentially 100% of the time even if the energy eigenstates 
that are being used in the experiment contain finite am- 
plitudes of the 'wrong' state [In particular, here we have 
in mind the states in Eq. ([32)) ]. We should emphasize 
here that present-day readout techniques are far from 
this ideal limit . 



B. Multi-qubit systems 

We now generalize some of the results concerning pos- 
sible errors in a typical experiment to the case where the 
system contains more than two qubits. The discussion 
below also demonstrates an interesting issue related to 
the significance of setting J — in the off state (here we 
have in mind the parametric coupling scheme). 

In order to clearly identify whether a given effect is 
related to mediated coupling or not, we start by consid- 
ering a multi-qubit system with no couplers (The expres- 
sions obtained below for the errors also apply if we inte- 
grate out the couplers in the circuit and use the effective, 
i.e. direct-plus-mediated, coupling strengths between the 
qubits). We first take the two-qubit Hamiltonian: 

A + + ^86™. (40) 

The above Hamiltonian has the interesting property that, 
with the proper assignment of the labels 00, 01, 10 and 
11 to the eigenstates of the Hamiltonian, the relations 

Eio — Eqq — En — Eqi 
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Eqi — Eqq — En — Eiq 
(Viol^^oo) = (^iil^^oi) 
(Voil^Voo) = (V'nl^lVio) 



(41) 



hold regardless of the values of the different parameters in 
the Hamiltonian. The above relations imply that the first 
two sources of error discussed in the previous subsection 
vanish completely for this system. It is therefore not 
required to have J12 < |Ai — A2I in order to perform an 
almost error-free quantum algorithm in this system. 

The question now is whether this situation persists for 
a system with more than two qubits. We approach this 
question by performing numerical simulations of a one- 
dimensional chain of three to ten qubits with the Hamil- 
tonian: 



N-l 

3=1 3 = 1 

From the numerical calculations we find that the gener- 
alized version of the relations in Eq. (j41j) holds only for 
the two qubits at the ends of the chain. The generalized 
version of the resonance-frequency relation [i.e. the first 
two lines in Eq. (|4"T]) ] continues to hold for all the qubits 
(i.e. the resonance frequency for flipping qubit j is inde- 
pendent of the states of all the other qubits). The Rabi 
frequency of qubit j oscillations, however, now generally 
depends on the states of the other qubits. 

In order to get an estimate for the above-mentioned 
errors, we take the three-qubit case and calcu- 
late the standard deviation 5f2R a bi.2 m the quantity 
~ (2) 

\&x \ipnOm)- When the coupling strengths Jjj+i 
are small compared to the differences between qubit gaps, 
we find that 



J 2 J 



23 



|(A 1 -A 2 )(A 2 -A 3 )(A 1 -A 3 ) 2 



(43) 



with a prefactor of order one (we also find that this ex- 
pression remains valid when we add a few more qubits 
at the ends of the chain). It is for this reason that one 
would like to set J = (for all pairs of qubits) in the off 
state, even though violating this requirement does not 
necessarily have detrimental effects on a two-qubit ex- 
periment |42|. It will also be highly desirable to set the 
gaps Aj to different values, even if they are not nearest 
neighbours. 

We now take a three-qubit chain with couplers, and we 
perform numerical simulations for several sets of param- 
eters with the leading-order coupling strengths (i.e. the 
parameters that correspond to J in the two-qubit case) 
set to zero. We calculate the resonance frequencies for 
qubit j (j = 1,2,3) for the four different states of the 
two other qubits. We find that the resonance-frequency 
shifts are consistent with the residual coupling strength 
given in Eq. (j2"5|) . We therefore expect our results con- 
cerning the residual coupling to hold for systems with 
more than two qubits. We also expect the form of 



the energy eigenstates to follow the rather straightfor- 
ward generalization of Eq. (|3"2"|) . If a distant pair of 
qubits have the same value of Aj, the energy scale as- 
sociated with any possible hybridization of energy levels 
(e.g. involving states of the form |00n20n30r^010r£6-") 
and 1 10n20n30n4000n6...)) will be small enough that it 
can be neglected. Similarly, hybridization between en- 
ergy levels that require flipping the states of several 
qubits should be greatly suppressed. 

Next we consider a loop of three or four qubits, i.e. in 
triangle and square geometries. Denoting the typical 
scale of the coupling strengths between neighbouring 
qubits by J and the detuning between the qubits by Aw, 
which is taken to be much smaller than the qubit fre- 
quencies, we find that the standard deviation in the res- 
onance frequency <5w rC sonance and the standard deviation 
in the transition matrix element <5f2R a bi when attempting 
to change the state of a given qubit are given by 



resonance 



<5^Rabi ~ 

in the three qubit case and 
5lu 1 

<5^Rabi 



Alu 2 

3 



^resonance 



J 

Alu 



J 

" Alu 3 
Alu 



(44) 



(45) 



in the four qubit case. The detailed expressions for the 
above quantities depend rather nontrivially on the values 
of the different parameters in the problem. Errors stem- 
ming from the spreads in resonance and Rabi frequencies 
therefore depend on the geometry of the multi-qubit sys- 
tem, e.g. chain versus closed loop. A two-dimensional 
lattice would belong to the latter category. 

Finally, we consider the geometry where every qubit 
is coupled to every other qubit, which is relevant to the 
case of a single coupler mediating coupling between sev- 
eral qubits. We find that the scaling of errors in this 
case follows the triangle geometry discussed above, with 
Aw representing the typical spacing of frequencies in the 
system. In particular, adding qubits to the system at dis- 
tant frequencies has little effect on the errors associated 
with a given qubit. 

The main result of this subsection is the fact that 
certain errors that have not been relevant to past ex- 
periments (i.e. two-qubit experiments) will likely arise 
in future experiments. Although this result might seem 
somewhat discouraging, it is encouraging that the errors, 
quantified by the spreads in the resonance and Rabi fre- 
quencies, follow scaling laws that make them less serious 
than one might intuitively expect (e.g. <5f?R a bi scales as 
the third or fourth power of J j Alu, as opposed to being 
linearly proportional to J / Alu). 
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VI. CONCLUSION 

We have analyzed the problem of a high-excitation- 
energy quantum object mediating coupling between two 
qubits. After reviewing some known results concerning 
the leading-order term in the mediated coupling, we ob- 
tained expressions that describe the residual-coupling in 
the off state and the entanglement in the energy eigen- 
states of the system. We have argued that our approach 
analyzing the properties of the eigenvalues and eigen- 
states of the total Hamiltonian is the appropriate one to 
describe recent, and possibly future, experiments. We 
have also estimated the expected errors originating from 
the non-ideality of the off state in typical experimental 
situations. Our results should be helpful in designing fu- 
ture circuits that employ the mediated-coupling approach 
in order to achieve tunable coupling between qubits. In 
particular, our results suggest that with properly chosen 
design parameters, the residual coupling in the off state 
could be greatly reduced in future experiments. 

We have focused on the case of a two-qubit system, 
which is the relevant case for present-day experiments. 
However, we expect our results for the errors in a two- 
qubit system to apply to multi-qubit systems as well. As 
we have shown, other sources of error that are absent in 
a two-qubit system arise for systems with more than two 
qubits. It would therefore be interesting and important 
in the future to analyze in more detail the properties of 
large, many-qubit systems. 

Finally, we should mention that although a large part 
of our analysis was formulated in the language of su- 
perconducting systems, our results are quite general 
and should apply to other systems that employ similar 
mediated-coupling mechanisms [43j . 

Appendix: Quasi-degenerate perturbation 
theory 



In this appendix we present a perturbation-theory pro- 
cedure for the case where the separation between some 
energy levels is not large compared to the energy scale of 
the perturbation [37| . 

Let us take the Hamiltonian 

H = H + V, (46) 

where we want to treat V as a perturbation. We assume, 
however, that the energy scale of V is not necessarily 
small compared to the energy separations within a subset 
of n eigenstates of Hq. The energy scale of V is small 
compared to the energy separation between this subset 
of levels and all the levels outside it. 

We now want to find the approximate energy levels and 
eigenstates of H in the vicinity of these n original energy 
levels of Hq. We proceed by assuming that these closely 
spaced eigenstates can have a large amount of mixing 
among themselves because of the added perturbation V , 
but we assume that mixing with all other states will be 
small. In other words, we express (any one of) the eigen- 
states of interest as 

n oc 

l*i)=£/«l^>+ E ( 47 ) 

3=1 j=n+l 

where \ipj) are the eigenstates of Hq with 

Hol^) = (48) 

the states \^>i) with j = 1, 2, n represent the states of 
interest, and gij are understood to be small enough to be 
treated perturbatively (this is the only reason we express 
the amplitudes /y and using two different symbols). 
We can now express the eigenvalue problem as 



oc I n 



£/yejM'j>+ E fe^l^l + E E 9i3VkMk)\=E i \Y,fii\^)+ E • (49) 



3=1 



j=n+l 



fc=l \ j=l 



j=n+l 



\j=i 



j=n+l 



If we multiply the above equation from the left by (tpi\ which can be rewritten as 



with 1 < I < n, we obtain the equation 

n oc 

fuei+J2fv V H+ E BijVij = Ekfu. (50) 

3=1 



9ij - 7; - 

&i — (-1 



Equation (|52p is now treated using a perturbative (or 
rather iterative) approach. If we neglect the sum over 
If we multiply the equation by (ipi\ with I > n, we obtain the states with j = n+ 1, we obtain: 



instead: 



YTj=i fij V ij 



gaei + E fi3 V l3 + zJ SijVij = Ei ga (51) 

3=1 j=n+l 



(53) 



Ei - ei 

Substituting this expression into Eq. (|50p , we obtain (for 
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every I with 1 < I < n): If we want to obtain more accurate results, we can 

take Eq. ([53]) and substitute it into the right-hand side 
« / 00 v V \ of Eq . (|52 |) . Expressing gj~ in Eq. (|53|) as ofi cv . we obtain: 

E ^ + V » + E 7— "p" & = /; ' / - ( 54 ) 

j=l y j''=n+l 1 J / 

where <5y is the Kronecker delta function. 



S"=i Aj^j 1 

ft; R 



— ej Ei - ei 



E *T% 



00 v^ 11 



' E Sr— ^ • (55) 



— e; -Ei -£(1. , J3j — £j 



Using this expression for gij in Eq. (|50|) . we find that to the next order: 



j=l V i'=n+l J j'=n+l j"=n+l y 3 J\ J 



The generalization to all orders is now obvious, if needed. 
One must be careful, of course, that the denominators 
on the left-hand side of Eqs. (f54|) and (f56|) contain the 
eigenvalue . The solution is therefore not yet obtained 
by a straightforward diagonalization of an n x n matrix. 
However, it is usually a good first approximation to use 
some (averaged) value for Ei on the left-hand side of the 
equation. This value can be taken from the energy levels 
of the original Hamiltonian Hq. One can then obtain 
more accurate results by introducing a second iterative 
procedure. Every time we obtain an approximate value of 
the energy Ei, we can substitute it into the left-hand side 
of Eq. ([5^| or to obtain an even more accurate result. 
Note that this iterative procedure is independent of the 
other one introduced above, i.e. the iterative substitution 
of Eq. ([53]) into Eq. (|5"!Z)) . In order to reach a given level 
of accuracy, both procedures must be performed to the 



appropriate order. 



ACKNOWLEDGMENTS 

We would like to thank M. Grajcar, J. R. Johansson, 
A. Maassen van den Brink, D. Tsomokos and J. Q. You 
for useful discussions. This work was supported in part 
by the National Security Agency (NSA), the Laboratory 
for Physical Sciences (LPS), the Army Research Office 
(ARO), the National Science Foundation (NSF) grant 
No. EIA-0130383 and the Japan Society for the Promo- 
tion of Science Core- To-Core (JSPS CTC) program. One 
of us (SA.) was supported by the Japan Society for the 
Promotion of Science (JSPS). 



[1] For recent reviews on the subject, see e.g. J. Q. You and 
F. Nori, Phys. Today 58 (11), 42 (2005); G. Wendin and 
V. Shumeiko, in Handbook of Theoretical and Computa- 
tional Nanotechnology, ed. M. Rieth and W. Schommers 
(ASP, Los Angeles, 2006). 

[2] Yu. A. Pashkin, T. Yamamoto, O. Astafiev, Y. Naka- 
mura, D. V. Averin, and J. S. Tsai, Nature 421, 823 
(2003). 

[3] P. R. Johnson, F. W. Strauch, A. J. Dragt, R. C. Ramos, 
C. J. Lobb, J. R. Anderson, and F. C. Wellstood, Phys. 
Rev. B 67, 020509(R) (2003). 

[4] A. J. Berkley, H. Xu, R. C. Ramos, M. A. Gubrud, F. W. 



Strauch, P. R. Johnson, J. R. Anderson, A. J. Dragt, C. 

J. Lobb, and F. C. Wellstood, Science 300, 1548 (2003). 
[5] T. Yamamoto, Yu. A. Pashkin, O. Astafiev, Y. Naka- 

mura, and J. S. Tsai, Nature 425, 941 (2003). 
[6] A. Izmalkov, M. Grajcar, E. Il'ichev, Th. Wagner, H.- 

G. Meyer, A. Yu. Smirnov, M. H. S. Amin, A. Maassen 

van den Brink, and A. M. Zagoskin, Phys. Rev. Lett. 93, 

037003 (2004). 

[7] H. Xu, F. W. Strauch, S. K. Dutta, P. R. Johnson, R. 
C. Ramos, A. J. Berkley, H. Paik, J. R. Anderson, A. J. 
Dragt, C. J. Lobb, and F. C. Wellstood, Phys. Rev. Lett. 
94, 027003 (2005). 



14 



[8] R. McDermott, R. W. Simmonds, M. Steffen, K. B. 

Cooper, K. Cicak, K. D. Osborn, S. Oh, D. P. Pappas, 

and J. Martinis, Science 307, 1299 (2005). 
[9] J. B. Majer, F. G. Paauw, A. C. J. ter Haar, C. J. P. M. 

Harmans, and J. E. Mooij, Phys. Rev. Lett. 94, 090501 

(2005). 

[10] B. L. T. Plourde, T. L. Robertson, P. A. Reichardt, T. 

Hime, S. Linzen, C.-E. Wu, and J. Clarke, Phys. Rev. B 

72, 060506(R) (2005). 
[11] M. Grajcar, A. Izmalkov, S. H. W. van der Ploeg, S. 

Linzen, T. Plecenik, Th. Wagner, U. Hiibner, E. Il'ichev, 

H.-G. Meyer, A. Yu. Smirnov, P. J. Love, A. Maassen 

van den Brink, M. H. S. Amin, S. Uchaikin, and A. M. 

Zagoskin, Phys. Rev. Lett. 96, 047006 (2006). 
[12] M. Steffen, M. Ansmann, R. C. Bialczak, N. Katz, E. 

Lucero, R. McDermott, M. Neeley, E. M. Weig, A. N. 

Cleland, and J. M. Martinis, Science, 313, 1423 (2006). 
[13] S. H. W. van der Ploeg, A. Izmalkov, A. Maassen van 

den Brink, U. Huebner, M. Grajcar, E. Il'ichev, H.-G. 

Meyer, and A. M. Zagoskin, Phys. Rev. Lett. 98, 057004 

(2007). 

[14] R. Harris, A. J. Berkley, M. W. Johnson, P. Bunyk, S. 

Govorkov, M. C. Thorn, S. Uchaikin, A. B. Wilson, J. 

Chung, E. Holtham, J. D. Biamonte, A. Yu. Smirnov, 

M. H. S. Amin, and A. Maassen van den Brink, Phys. 

Rev. Lett. 98, 177001 (2007). 
[15] T. Hime, P. A. Reichardt, B. L. T. Plourde, T. L. Robert- 
son, C.-E. Wu, A. V. Ustinov, and J. Clarke , Science 

314, 1427 (2006). 
[16] A. O. Niskanen, K. Harrabi, F. Yoshihara, Y. Nakamura, 

and J. S. Tsai, Phys. Rev. B 74, 220503(R) (2006). 
[17] A. O. Niskanen, K. Harrabi, F. Yoshihara, Y. Nakamura, 

S. Lloyd, and J. S. Tsai, Science 316, 723 (2007). 
[18] J. Plantenberg, P. C. de Groot, C. J. P. M. Harmans, 

and J. E. Mooij, Nature 447, 836 (2007). 
[19] M. A. Sillanpaa, J. I. Park, and R. W. Simmonds, Nature 

449, 438 (2007). 
[20] J. Majer, J. M. Chow, J. M. Gambetta, J. Koch, B. R. 

Johnson, J. A. Schreier, L. Frunzio, D. I. Schuster, A. 

A. Houck, A. Wallraff, A. Blais, M. H. Devoret, S. M. 

Girvin, R. J. Schoelkopf, Nature 449, 443 (2007). 
[21] D. V. Averin and C. Bruder, Phys. Rev. Lett. 91, 057003 

(2003) . 

[22] B. L. T. Plourde, J. Zhang, K. B. Whaley, F. K. Wilhelm, 
T. L. Robertson, T. Hime, S. Linzen, P. A. Reichardt, 
C.-E. Wu, and J. Clarke, Phys. Rev. B 70, 140501 (R) 

(2004) . 

[23] M. Wallquist, J. Lantz, V. S. Shumeiko, and G. Wendin, 
New J. Phys. 7, 178 (2005). 

[24] A. Maassen van den Brink, A. J. Berkley, and M. 
Yalowsky, New J. Phvs. 7, 230 (2005). 

[25] Before Refs. [2l|, [22l, [24| there had been proposals to cou- 
ple superconducting charge qubits using a large induc- 
tance; see, e.g., A. Shnirman, G. Schon, and Z. Hermon, 
Phys. Rev. Lett. 79, 2371 (1997); Y. Makhlin, G. Schon, 
and A. Shnirman, Nature 398, 305 (1999). Smaller, more 
realistic, inductances were considered in J. Q. You, J. S. 
Tsai, and F. Nori, Phys. Rev. Lett. 89, 197902 (2002). 
A similar approach for flux qubits was proposed in J. Q. 
You, Y. Nakamura, and F. Nori, Phys. Rev. B 71, 024532 

(2005) . 

[26] See also M. D. Kim and J. Hong, Phys. Rev. B 70, 184525 

(2004); M. D. Kim, Phys. Rev. B 74, 184501 (2006). 
[27] E. Farhi and S. Gutmann, Phys. Rev. A 57, 2403 (1998). 



[28] P. Bertet, C. J. P. M. Harmans, and J. E. Mooij, Phys. 

Rev. B. 73, 064512 (2006). 
[29] A. O. Niskanen, Y. Nakamura, and J. S. Tsai, Phys. Rev. 

B 73, 094506 (2006). 
[30] M. Grajcar, Y.-X. Liu, F. Nori, and A. M. Zagoskin, 

Phys. Rev. B 74, 172505 (2006). 
[31] Performing two-qubit gates on superconducting qubits 

by driving the qubits at the sum or difference of their 

characteristic frequencies was first proposed in Y.-X. Liu, 

L.-F. Wei, J. S. Tsai, and F. Nori, Phys. Rev. Lett. 96, 

067003 (2006). 

[32] Experiments on a single qubit coupled to a 'cavity' were 
first reported in I. Chiorescu, P. Bertet, K. Semba, Y. 
Nakamura, C. J. P. M. Harmans, and J. E. Mooij, Nature 
431, 159 (2004); A. Wallraff, D. I. Schuster, A. Blais, L. 
Frunzio, R. S. Huang, J. Majer, S. Kumar, S. M. Girvin, 
and R. J. Schoelkopf, Nature 431, 162 (2004); See also J. 
Johansson, S. Saito, T. Meno, H. Nakano, M. Ueda, K. 
Semba, and H. Takayanagi, Phys. Rev. Lett. 96, 127006 

(2006) . 

[33] There are a large number of theoretical studies on this 
subject. See, e.g., J. Q. You and F. Nori, Phys. Rev. B 
68, 064509 (2003); R. Migliore, A. Konstadopoulou, A. 
Vourdas, T. P. Spiller, and A. Messina, Phys. Lett. A 
319, 67 (2003); Y.-X. Liu, L. F. Wei, J. R. Johansson, 
J. S. Tsai, and F. Nori, Phys. Rev. B 76, 144518 (2007); 
A. Blais, J. Gambetta, A. Wallraff, D. I. Schuster, S. 
M. Girvin, M. H. Devoret, and R. J. Schoelkopf, Phys. 
Rev. A 75, 032329 (2007); H. Nakano, K. Kakuyanagi, 
M. Ueda, and K. Semba, Appl. Phys. Lett. 91, 032501 

(2007) ; F. Helmer, M. Mariantoni, A. G. Fowler, J. von 
Delft, E. Solano, F. Marquardt. iarXiv:0706.3625l 

[34] C. Hutter, A. Shnirman, Y. Makhlin, and G. Schon, Eu- 
rophys. Lett. 74, 1088 (2006). 

[35] See, e.g., G. Baym, Lectures on Quantum Mechanics 
(Addison- Wesley, New York, 1990). 

[36] We shall not derive the expressions for Cj in this section; 
the relevant expressions agree, to lowest order, with those 
that will be derived in Sec. IV. 

[37] See, e.g., J. J. Sakurai, Modern Quantum Mechanics 
(Addison- Wesley, New York, 1994). 

[38] In the theoretical treatment of this system, one can also 
rewrite the Hamiltonian such that the matrix elements 
Aoo and -Boo are absorbed into ei and £2- One then re- 
places An by An — Aqo, and similarly for B. 

[39] The second term on the right-hand side of Eq. (1321) is 
obtained by including the gaps Aj in the denominator 
in Eq. (|20|) . which slightly modifies Ji depending on the 
state under consideration. 

[40] M. Wallquist, V. S. Shumeiko, and G. Wendin, Phys. 
Rev. B 74, 224506 (2006). Note that, strictly speaking, 
the coupler is no longer an harmonic, i.e. linear, oscillator 
in this design. 

[41] For example, it is possible, in principle, for a device 
that couples to the operator a z of a two-level system 
to distinguish with essentially 100% fidelity between 
energy eigenstates of the form \/0.6|0) + V 0.4| 1) and 

Voa\o) - VoI\i). 

[42] For coupling proposals that rely on the relation J < 
Ai — A2I in order to achieve effective decoupling, see 
e.g. C. Rigetti, A. Blais, and M. Devoret, Phys. Rev. 
Lett. 94, 240502 (2005); G. S. Paraoanu, Phys. Rev. B 
74, 140504(R) (2006); S. Ashhab, S. Matsuo, N. Hatak- 
enaka, and F. Nori, Phys. Rev. B 74, 184504 (2006); S. 



15 



Ashhab and F. Nori, Phys. Rev. B 76, 132513 (2007); 
See also Ref. [3lj . 
] For example, there has been a recent proposal to use 
the high-frequency electrons as mediators for performing 



gates on nuclear spins in NMR quantum computing; J. 
S. Hodges, J. C. Yang, C. Ramanathan, and D. G. Cory, 
arXiv:0707.2956vl. 



